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The recently formulated model of highly-anisotropic and strongly dis- 
sipativc hydrodynamics is used in 3+1 dimensions to describe flow char- 
acteristics and strangeness production in Au+Au collisions at the highest 
RHIC beam energy. Our results show very weak dependence on the initial 
momentum anisotropy, provided the anisotropic phase lasts no longer than 
1 fm/c. 

PACS numbers: 25.75.-q, 25.75.Ld, 24.10.Nz, 25.75.Nq 

1. Introduction 

Soft-hadronic data collected in the ultra-relativistic heavy-ion experi- 
ments may be well described in the framework of the standard perfect-fluid 
hydrodynamics pQ or by dissipative hydrodynamics with small viscosity 
[21 El HI El El [7] . However, the success of those approaches relies strongly 
on the assumption that the produced system reaches state of local thermal 
equilibrium within a fraction of a fermtU. 

Most of the microscopic models of the early stages fail to explain such 
short thermalization times pQ. This difficulty is known as the early ther- 
malization puzzle. One of its common solutions is the concept of a strongly 
coupled quark-gluon plasma [8]. In addition, many microscopic approaches 
assume that the produced system exhibits initially large anisotropies in the 
momentum space, for example, see [9]. 

Recently, several models have been developed [HUE], which include the 
early highly-anisotropic phase of the collisions. This has been achieved intro- 
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ducing a pre-equilibrium stage connected with a subsequent perfect-fluid de- 
scription. Combination of different approaches in a single framework seems, 
however, not completely satisfactory. The need for a concise model which 
can describe different stages of heavy-ion collisions in the uniform way trig- 
gered our development of the highly-Anisotropic and strongly-Dissipative 
HYDROdynamics (ADHYDRO) [13 EJ O [E] , see also pH [T7] 

In this paper we present our results obtained within the ADHYDRO 
framework coupled to THERMINATOR [TBI 02]. For the first time we use 
ADHYDRO in 3+1 dimensions. We discuss flow characteristics of the emit- 
ted hadrons and strangeness production depending on the initial momentum 
anisotropy. 

2. Formulation of the model 

In the ADHYDRO model the evolution of the system is described by 
the following equations [12] 

d^ v = 0, (1) 
V = S ' ( 2 ) 

which express the energy-momentum conservation and the entropy produc- 
tion laws, respectively. The energy-momentum tensor T^ u in Eq. ([1]) has 
the form 

T" u = (e + P ± ) U^U V - P ± - (P ± - P\\)V X V U , (3) 

which allows for the asymmetry between the longitudinal, Pm, and trans- 
verse, P±, pressures. In the limit where the system becomes isotropic, 
Pi I = Pj_ = P, the formula reproduces the energy-momentum tensor 
of the perfect fluid. Similarly, the entropy production law ([2]) is reduced 
to the entropy conservation law, if we assume that the entropy source, S, 
vanishes. The four-vector defines the four-velocity of the fluid and 
is the four-vector defining the beam axis. In the general case, and 
may be parametrized in the following way 

W i = (uq coshi?, u x , u y , uo sinh#), (4) 
V = (sinh 0, 0,0, cosh •&), (5) 

where u x and u y are the transverse components of the four-velocity veloc- 
ity field (u± = \Ju"i + u 2 and uq = <Jl + u\), $ is the longitudinal fluid 

rapidity. The parametrizations @ and ([5]) satisfy simple normalization con- 
ditions, i.e., U 2 = 1, V 2 = -1, U ■ V = 0. The entropy flux a» in © is 
defined by the formula 

a fl = aU f ", (6) 



sqm2011 printed on November 9, 2011 



3 



where a is the non-equilibrium entropy density. 

One can show [12] that instead of P» and P± it is more convenient to use 
the entropy density a and the anisotropy parameter x as two independent 
variables (to a good approximation we have Pn/P± = x~ 3 / 4 ). Similarly 
to standard hydrodynamics with vanishing baryon chemical potential, the 
energy density e introduced in Eq. ([3|), the entropy density a, and the 
anisotropy parameter x are related through the generalized equation of state 
£ = e(a,x). In our model we use the following ansatz [J3] 

s(x,a) = e qgp (o-)r(x), (7) 
P±(x,a) = P qgp (a) [r(x) + 3xr'(x)} , 
P\\(x,a) = iqgp(cr) [r(x) - 6xr'(x)] . 

where £ qgp and P qgp define the realistic equation of state constructed in Ref. 
[20] . The function r(x) is the pressure relaxation function characterizing the 
properties of the fluid which exhibits the anisotropy x. Here we use the form 
introduced in fT2l 



r(x) 



1 



x 



arctan \/x — 1 



(8) 



In the isotropic case x = 1, r(l) = 1, r'(l) = 0, and Eq. ([7]) is reduced to 
the equation of state used in 



The function E which appears on the right-hand-side of Eq. ([2]) defines 
the entropy production due to microscopic processes taking place in the 
system. Exactly these processes lead to thermalization of the system. We 
use the form proposed in p2] 

E(c7,x) = -= , (9) 



where the time-scale parameter r eq = 0.25 fm controls the rate of equili- 
bratior|l|. In the limit of small anisotropy Eq. ([9]) is consistent with the 
quadratic form of the entropy production in the Israel-Stewart theory. Far 
from equilibrium, hints for the form of E are lacking, although we may 
expect some suggestions from the AdS/CFT correspondence [21] . Thus, 
for large anisotropies the formula ([9]) should be treated as an assumption 
defining the dynamics of the system. 



2 Using this value we find that the system equilibrates within about 1 fm. 
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3. Initial conditions and freeze-out 

In the general 3+1 case we have to solve Eqs. JT]) and (J2J for five 
unknown functions a, x, u x , u y , and v), which depend on the space-time 
coordinates: r, xj_, and r/ (r is the proper time and rj is the space-time 
rapidity). Since the system's evolution is treated hydro dynamically from the 
very early stages where the anisotropics are expected to be very large, we fix 
the initial starting time for ADHYDRO to tq = 0.25 fm. Similarly to other 
hydrodynamic calculations, we assume that there is no initial transverse 
flow, u x (to , xj_ , rj) = u v (to, xj_, 77) = 0. For the initial longitudinal rapidity 
of the fluid we assume the Bjorken scaling $(to, xj_, rj) = rj. We check 
three scenarios: i) the initial source is strongly oblate in the momentum 
space, x(tq, xj_, rf) = 100, which corresponds to a transversally thermalized 
source, ii) the source is prolate in momentum space, x(tq, xj_, rf) = 0.032 
(this case is analogous to i) because r(100) = r(0.032)), and iii) the initial 
source is isotropic, x(to, xj_, 77) = 1, which gives the closest description to 
the standard hydrodynamic case. The initial entropy density profile is given 
by the formula 



o"(t-o,x_l,?7) = e gq p 



—— rr , p(b,X ± ,T]) = , 10 

r{x{T ,x ± ,rj))\ p(0,0,0) 

where p is the normalized initial density of sources. The density profile p is 
given as the tilted source worked out in Ref. [22J and e~^ p is the inverse of 
the function e gq p(a). The initial energy density £\ = 107.5 GeV/fm 3 is the 
same for all three analyzed cases. 

The evolution is determined by the hydrodynamic equations until the 
entropy density drops to o"f = 1.79 fm -3 , which for x = 1 corresponds to 
the temperature Tf = 150 MeV. According to the single-freeze-out scenario, 
at this moment the abundances and momenta of particles are expected 
to freeze-out and particles freely stream to detectors. The processes of 
particle production and decays of unstable resonances are described by using 
THERMINAT0R 2 [19], which applies the Cooper- Frye formalism to generate 
hadrons on the freeze-out hypersurface extracted from ADHYDRO. 



4. Results and conclusions 

The THERMINATDR 2 package [19] allows us the calculation of several one- 
and two-particle observables (such as the particle pt spectra, the directed 
and elliptic flow coefficients, the HBT radii, etc.). Due to limited space, in 
this paper we show only the hyperon spectra and the directed flow coefficient 
V\. More results will be presented and discussed in a separate publication. 

The scaling (fTUj) helps us to keep the final particle multiplicities approx- 
imately the same in the three considered cases. In particular, this can be 
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Fig. 1. Transverse- momentum spectra of hiperons (left part) compared to the ex- 
perimental data from STAR [23J and directed flow coefficient (right part) compared 
to the STAR data [2J]. Theoretical lines are obtained from AD HYDRO for Au+Au 
collisions at ^/snn = 200 GeV and the centrality class c = — 5%. The results 
are shown for three cases: oblate source (blue dashed line), prolate source (red 
dashed-dotted line) and perfect fluid (green dotted line) as discussed in Sect. [3J 



concluded from the left-hand-side of Fig. [T] where we show the transverse 
momentum spectra of hyperons. We observe that the stronger transverse 
flow in the case i) results in a bit harder spectra as compared to the case 
ii). Despite this fact, the spectra do not differ significantly. We observe 
that the model spectra of S's and fTs agree well with the data, while the 
normalization of A's is too small. 

The authors of Ref. [25] have shown that the v\ coefficient may be treated 
as a probe for measuring the thermalization time, since it is very sensitive to 
the early difference of pressures. On the right-hand-side of Fig. [I] we present 
our results for the directed flow coefficient v\. We observe small sensitivity 
to large initial anisotropy in the midrapidity region provided the anisotropic 
stage lasts not longer than 1 fm. Our results are not so much restrictive as 
the results presented in [25], since we do not fix the final multiplicities but 
allow them to vary within the experimental errors. 

More results obtained in the ADHYDRO model in 3+1 dimensions will 
be presented in a separate publication. We note that our previous results 
obtained in the 2+1 version [15J (a boost-invariant version) show similar, 
weak dependence of other physical observables on the initial pressure asym- 
metry. In conclusion, we state that we have found further evidence that 
there is a place for a highly-anisotropic phase at the early stages of heavy- 
ion collisions, provided the system reaches local thermal equilibrium before 
1 fm. 
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